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Jacobi elliptic functions are flexible functions that appear in a variety of problems in physics and 
engineering. We introduce and describe important features of these functions and present a physical 
example from classical mechanics where they appear: a bead on a spinning hoop. We determine 
the complete analytical solution for the motion of a bead on the driven hoop for arbitrary initial 
conditions and parameter values. 
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II. INTRODUCING JACOBI ELLIPTIC 
FUNCTIONS 
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In 1788, James Watt introduced a mechanical device 
into steam engines that moderated the flow of steam 
and controlled the speed of the engine^ This device, 
known as the governor, is composed of two solid pendula 
(weights) fixed to a common vertical axis. The power of 
the steam engine rotates the axis and lifts the weights 
to a height determined by the angular speed. The be- 
havior of the weights is nowadays presented to students 
as an example of Lagrangian mechanics as a mass con- 
strained to move frictionlcssly on the surface of a sphere 
with constant azimuthal frequency or a bead moving on 
a rotating hoop4^ 

Solving the equation of motion of the bead on a ro- 
tating hoop is generally avoided as it does not have a 
straightforward solution. The problem is usually limited 
to the determination of the stationary position of the 
mass for a given angular frequency of the hoop, although 
some papers treat the complete problem numerically or 
within some approximations.— ~— Here, we present the full 
analytical solution for the motion of the bead on the hoop 
which requires the use of Jacobi elliptic functions (JEFs). 

There are a number of physics problems where JEFs 
provide an excellent approximation if not a complete de- 
scription of the system.^ These special functions are 
also amply used in the conformal mapping of engineer- 
ing problems.— ^ JEFs have some interesting features. 
For example, trigonometric and hyperbolic functions are 
special cases of Jacobi elliptic functions. If considered in 
the complex plane, JEFs are doubly periodic, a feature 
appearing, for example, in bifurcation theory of chaotic 
systems 



Despite the diversity of physical problems where JEFs 
appear, these functions arc generally not part of the 
physics curriculum. This motivates us to offer a con- 
cise and pedagogical exposition of these functions. We 
then demonstrate how these special functions appear in 
the example of the bead on the hoop. 



In this section, we introduce Jacobi elliptic functions 
through elliptic geometry and integral inversion,^^— 
and discuss the geometrical interpretation of their argu- 
ments. 



Definitions and relations between Jacobi elliptic 
functions 



We start by redefining the basic trigonometric func- 
tions sine and cosine in terms of the functional inverse 
of specific integrals. Usually, these functions are intro- 
duced using circular geometry. That is, for a given point 
(x, y) on a circle of radius r, sin = y/r and cos = x/r. 
An alternative way to define these functions is to real- 
ize that each inverse trigonometric function is a solution 
of a definite integral. For example, simple trigonometric 
substitutions lead to 



dt 

VT - F 



= arcsmy, 



dt 



1 + t 2 



arctany, (1) 



with the conditions that — 1 < y < 1 and yl — t 2 > 0. 
We can reinterpret these relations by stating that they 
define the inverse trigonometric functions. Following 
Abel we obtain trigonometric functions by inverting these 
integrals For example, defining the argument 6 as 



0(y) = 



dt 



arcsmy, 



(2) 



we obtain the trigonometric sine function sin 6 = y. With 
the Pythagorean relation, we can define the cosine func- 
tion and all others follow. This alternative definition im- 
plies that the integrals 9{y) must in particular display the 
parity and periodicity of the corresponding trigonometric 
functions y{9). 

Consider now the case of elliptic geometry. Gauss, 
Legendre, Abel, Jacobi and Weierstrass studied in depth 
properties of elliptic functions and integrals. For the pur- 
pose of this paper we note only the result, proved by 
Legendre^ that any integrated expression containing a 
third or fourth degree polynomial in the denominator of 
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a fraction can be reduced to a linear combination of the 
following elliptic integrals of the first, second, and third 
kind, 



dt 




o v/a-^Hi-fc 2 * 2 ) 

o \f \ - k 2 sin 2 4>' 
dt 
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E((f>,k)= I Jl-k 2 sin 2 <j)' d<p' 
Jo 

U(y,k,n) 



dt 



o (1 - n 2 t 2 ) y/(l - t 2 ) (1 - k 2 t 2 ) 



II ((/), k, n) = 



o (l - n 2 sin 2 (f>') y/l - k 2 sin 2 </>' 



(4) 



(5) 



where y = sine/), i = sm<f>', k G (—1, 1) is the modulus 
and neR the characteristic. We assume all square roots 
to be positive. In Eqs. ([3])-([5]), the first expression (in t) 
is known as the Jacobi form whereas the second (in <f>') 
is Legendre's form. When = 7r/2, written F (■§,&) = 
K (fc), the integrals are said to be complete and otherwise 
incomplete. We will use the first two elliptic integrals 
exclusively in what follows. 

Because trigonometric functions are sometimes more 
useful than their inverses, it may be worthwhile to in- 
vert the elliptic integrals as well.— ^ In this process, Abel 
and Jacobi followed the work of Legendre and introduced 
what are now known as the Jacobi elliptic functions. 
These functions result from the inversion of the elliptic 
integral of the first kind. Introducing the Jacobi nota- 
tion u(4>) = F((j>, k) for this integral, we can formally 
invert the Legendre form of the integral, Eq. ((3]), to ob- 
tain 4>(u) = amu, the amplitude of u. In the Jacobi form, 
u(y) = F(y,k), the inversion of Eq. ((3]) leads one to de- 
fine the sine-amplitude Jacobi elliptic function written in 
Glaisher's notation (y <G [-1, 1], <f) € [— 7r/2, 7r/2]) 



sn(it, k) = y = sin<^ = sin (am it) 



(6) 



Note that the fc-dependence is often implicit in the liter- 
ature, resulting in the notation snti instead of sn(u, k). 
The last two equalities in Eq. ([6|) relate the Legendre and 
Jacobi forms of the elliptic integrals and justify the name 
given to this Jacobi elliptic function. Because Eq. ((3j) is 
odd, u(y) — F(y,k) = sn~ 1 (y,k) is odd in y as well. 
That is, the sine-amplitude is an odd function of u. 

It is natural to introduce two further Jacobi elliptic 
functions. Taking the cosine of the amplitude of it, we 
define the cosine-amplitude 




K(k) 2K(k) 3K(k) 4K(k) 



FIG. 1. (Color online) Representative graphs of the Jacobi 
elliptic functions sn(w), cn(u), and dn(u) at fixed value of the 
modulus k = 0.9. In the bead and the hoop problem, that 
argument is a function of time. 



Finally, in the Legendre form all elliptic integrals in 
Eqs. ©-© contain A(<f>) = y/l - k 2 sin</>. This leads 
one to define the delta-amplitude JEF that can also be 
written as a derivative of the amplitude^ 



dn(w, k) — vT— k 2 sn 2 u 



d(axrm) 
du 



(8) 



Equation can also be used to obtain the last equality. 
Figure [1] displays the characteristic behavior of the sine-, 
cosine- and delta-amplitude functions. 

From the definitions above immediately follow the 
identities 



sn 2 u + cn 2 u 



1. 



dn it + k 2 sn 2 u = 1, 
cn 2 u + (1 — k 2 ) sn 2 u = dn 2 tt. 



(9) 
(10) 
(11) 



cn(w, k) = cose/) = cos (amu) 



(7) 



For limiting values of fc, the Jacobi elliptic functions 
reduce to trigonometric and hyperbolic functions 

k = : snw = sin u, cnu = cos u, dnu — 1, (12) 
k = 1 : smt = tanh u, cnu = dnu = sechtt. (13) 

Jacobi elliptic functions thus include trigonometric and 
hyperbolic functions as special cases. However, JEFs are 
more than a simple generalization of elementary func- 
tions as can be seen when studying the properties of these 
functions in the complex plane (see Appendix [A"| . 

Finally, we mention that from the definitions 
one can define quotients of JEFs. One defines sc it = 
snu/cnit which carries the name tangent-amplitude 
scu = tn it. Other quotients are, for example, nd it = 
1/dn u, sdu — snu/dnit, etc. Further properties of JEFs 
that are used in the treatment of the bead on the hoop 
are summarized in Appendix [XJ 



B. Geometric interpretation of the arguments (it, k) 



An alternative way to define this function is to write 
cnu = \/l sn 2 u since, from Eq. ©, y € [—1,1]. 



In the case of trigonometric functions, such as y = 
sin 8, the argument 6 has an unequivocal geometric in- 
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tcrprctation in circular geometry; it is the angle between 
the x-axis and the radius at a specific point (x, y) on the 
circle. The question arises whether in a JEF the argu- 
ment u and the modulus k have a geometric meaning 
as well. We first caution that the argument of a special 
function (Bessel, Hypergeometric, etc.) does not gener- 
ally have such unequivocal physical or geometric mean- 
ing. For JEFs, the geometric definition of the argument 
k is unique but not for the argument u, an important 
point not discussed explicitly in the literature. 

The modulus k is given in elliptic geometry by k = e, 
the eccentricity of the ellipse e 2 = 1 — b 2 /a 2 , where a > b 
are the semi-axes (see Fig. [2} 



y_ 

b 2 



= l. 



(14) 



The equivalence between k and e is found by realizing 
that the arc length sbp on the ellipse is expressed in 
terms of the elliptic integral of the second kind. Using 
the parametrization x = a sin <f>, y = b cos <f> 1 the in- 
finitesimal arc length along the ellipse is 



ds 



dx 2 + dy 



a' cos* 



b 2 > 



= a\ 1 — e 2 sin 



(15) 



Using Eq. ((4]), the arc length sbp oi the ellipse from B 
to P can be written 



sbp = a 



1 -e 2 



h> A A' — 



aE(<j),k = e). (16) 



This relation between the arc length of the ellipse and 
the function E(<j>, k) is the reason why Fagnano (cited in 
Ref.H!) called F, E, and II elliptic integrals and provides 
a direct interpretation of the modulus of the JEFs. 

The geometric interpretation of JEFs' argument u is 
less direct and not unique. Fagnano discovered that for 
k^ 1 = a/2, the elliptic integral of the first kind describes 
the arc length on the lemniscate defined by the relation 
r 2 = cos 28 (using the same definitions of r and 8 as 
in Fig. H])] 12 ! 14 Serret, building on the work of Fagnano 
and Legendre, found a set of curves on which F(jf), k) at 
arbitrary value of k £ (—1,1) is an arc length (see, e.g. 
Ref . [l2l) . There is also a counterpart in three dimensions: 
Seiffert's spirals defined on a sphere. The arc length of 
the spiral is given by the first elliptic integral £ In these 
examples the argument u of the JEF is the length of the 
Serret curve or the Seiffert spiral. 

Additional geometric interpretations of the argument u 
can be found from trigonometry as it applies to an ellipse 
using geodesy and elliptic geometry. This is discussed in 
detail in Refs.l5landli~7l. All of these interpretations relate 
u to a length of some geometric form related to an ellipse. 
It is important, however, to realize that the angles used 
to define u are different for each physical situation. This 
leads to the important insight that the argument u has 
no unique geometric interpretation beyond being an arc 
length and only has a physical interpretation in some 
special cases^ i 14 ' 17 




FIG. 2. Representations of points and angles on the ellipse 
that appear in the definitions and relations involving Jacobi 
elliptic functions. The figure shows the polar angle 6 and 
the complementary of the eccentricity angle <j>. The smaller 
circle of radius 6 (dashed line) is inscribed in the ellipse and 
has common points at ±_B = (0, ±6, 0). The ellipse is also 
inscribed in the larger circle with radius a (dash-dotted line) 
which has two points in common ±A = (±a, 0, 0). 



III. THE BEAD ON THE HOOP 

In this section we present an application of Jacobi ellip- 
tic functions to a well-known problem. Consider a bead 
of mass m that slides without friction on a hoop of ra- 
dius R rotating with constant angular velocity w = <p (see 
Fig. IS]). In spherical coordinates (r,8,4>) the Lagrangian 
C = T — V for this system is 



£ 0, 8 ,t 



1 



uj 2 sin 2 1 



-Mi 



COS! 



, (17) 



where Iq = raR 2 and Vq = mgR. The zero point of the 
potential energy is located at the bottom of the hoop. 
The angles 8 and <j> should not be confused with the an- 
gles of same name in Sec. [TTJ The polar angle 8 is mea- 
sured from the bottom of the hoop for consistency with 
the convention usually used in the limiting case of the 
simple non-linear pendulum. 

There are three variables (r, 8, <f>) and two constraints 
r = R and <f> = cj, implying that there is only one gener- 
alized coordinate, 8. That is, the motion of the bead is 
described by the time-dependent function 8 (t) . 

The discussion of the physical content of the follow- 
ing calculations is greatly aided by the introduction of 
an effective potential V c s. Since 8 is the only degree of 
freedom, we combine the azimuthal kinetic energy with 
the gravitational potential energy and rewrite Eq. f| 1 T[) 

as C (6, 8, t] = T g (e) - V cff (8), with 



Vcs (g) 
V 



--l — \sm z 8 + (l-cos8) (18) 



FIG. 3. Schematic of a bead on a hoop rotating about the 
vertical z-axis at a constant azimuthal frequency uj = <f>. The 
diamond at 9 — ±9 S locates the global minimum of the ef- 
fective potential at the commonly named stationary point, 
whereas the star at 9 — ±9 e locates the so-called ejection 
point (see text and Fig. 3]). The positions of 9 S and 9 e are 
obtained using uj = 3oj c /2. The motion of the bead is de- 
scribed by the polar angle 6(t) measured from the bottom of 
the hoop. 



FIG. 4. (Color online) The effective potential (|18p as a func- 
tion of 9. The dashed line is for uj < uj c (for the present plot 
we chose u) = 0J c /5), the dotted line for uj — uj c , and the solid 
line is for uj > uj c (here uj — 3w c /2). The diamond at ±# s 
and star at ±9 e are indicated in Fig. [3] for the present set of 
parameters. The effective potential energy at 9 e and 9 = is 
always zero. Note that VSf* = V e s(n) = 2Vo- 



A. General Solution of the Euler-Lagrange 
Equation 



and uj 2 = Vq/Iq = g/R. As shown in Fig. [4j the qual- 
itative shape of the effective potential depends on the 
azimuthal frequency uj. One distinguishes two regimes: 
uj < uj c when 9 < 0, and uj > uj c when 9 > 0. Since 
9q = 7r is always a maximum of V e s, this position is un- 
stable. On the other hand, 9q = is a minimum in the 
first regime but a local maximum in the second regime 
where two new stable minima appear at ±(9 S . For en- 
ergies near the minima at ±9 S , the bead will oscillate 
about 9 S , remaining exclusively on one side of the hoop. 
The stable equilibrium point 9 S can be obtained by set- 
ting the derivative of the effective potential equal to zero. 
We find 



(19) 



which is valid in the supercritical regime (uj > uj c ). 

The differential equation governing 9(t) is obtained 
from the Euler-Lagrange equation associated with 
Eq. {3U 



— sin26>- ■ - 2 



cj: sin ( 



(20) 



and introduce the initial conditions 6(t = 0) = 9q and 
9 (t = 0) = 9q. Without loss of generality, we make use 
of the symmetry of the problem and define the initial an- 
gle 9q only on one half of the hoop #o G [0,7r]. Though 
nonlinear, this differential equation can be solved analyt- 
ically. The presence of sin 9 and sin 29 leads to solutions 
involving JEFs described in Sec. ILTl 



We begin by multiplying Eq. (|20[) by 9 and integrating 
to obtain 

-At- + cos 2 9(t) - 2cos9 s cos9(t) =C, (21) 



where the constant C is given by 



C = cos 2 O - 2cos6» s cos6» + f — j , (22) 



and 9 S is given in Eq. (|T9)) . Note that Eq. (|21j) is not 
valid for uj = 0; this situation must be considered sepa- 
rately (see Sec. lIIIB"3]) . It is worth mentioning that 9 S is 
nonzero only for uj > uj c and takes values in the interval 
9 S G [0, 7r/2]. For uj < uj c , only 6q — is a stable station- 
ary solution. It is then natural to extend the definition 
of 9 S by imposing 9 a = for the case where uj < uj c . 

Because C is a constant, the left hand side of Eq. (|2Tj) 
is a constant of motion of the system. The expression 
for C is of particular interest because it allows us to re- 
duce arbitrary initial conditions to ones with zero initial 

speed. Consider arbitrary initial conditions (9o,9q) for 
which 9 ^ and E tot < V"g ax . From Eq. ([H]) we can 
find a time to and new initial conditions (9q,9q) with 
zero initial angular speed 9o ~ such that 6 (to) = 6q 

and 9(to) = 9q. The new initial condition (#0j0) written 
in terms of the actual one (9q, 9q) is 



cos 9 — cos 9 S ± 



cos # — cos 9 S 



1 ~ 

UJ 



2 

2"0- 



(23) 
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Thus, arbitrary initial conditions with 9o / can be 
written in terms of initial conditions with (9 = 0. In the 
following, we limit our discussion to the motion of the 
bead on the hoop with zero initial speed unless E to t > 
7™. This case will be discussed later. 
We now use the Wcicrstrass substitutio: 



,18 



tan - 
2 



1 



COS ( 



1 + cos 



in Eq. (f2"T|) and obtain 



r(t) 



2dr' 



\[p~^ 



2(C + 1)t' 2 +p+ 



(24) 



(25) 



through separation of variables. This type of equation is 
usually studied numerically or in limiting cases only£~— 
Here we extend the calculation and invert the integral to 
obtain 9{t). The structure of the integral implies that the 
inverse is expressed in terms of JEFs defined in Sec. HU 

For convenience, we define To = tan ($o/2) and the two 
constants 



p± = C - 1 ± 2cos( 



(26) 



After factoring out p_ from the denominator of the in- 
tegrand in Eq. (|25[) . we write the remaining polynomial 
under the square root in the form 



r ,4 + 2 



C + l 



P+ 



P- J \P 
where the roots are given by 

7i 2 2 G W> b } 

with 



/2 2 
= I T - 7 X 



(27) 
(28) 



^ = 1 - cos go __ 1 + cos 9 - 2 cos S 

a ~ 1 + cos O ~ T ° ' ~ 1 - cos O + 2 cos S ' ^ ' 

The root a is always non-negative while b can take any 
real value. As a result, the relative values of the two 
roots a and b divides the range 0q G [0, tt] into one (for 
uj < lu c ) or two (for uj > uj c ) intervals. When uj > u c , 
the two intervals are divided by the point b = 0, which 
occurs when the initial condition 0q takes the value 



e (w) = arccos (2 cos S — 1) = arccos 2 



1 



(30) 



For 9o < 8 e we have b > while for 0q > 9 e we have 
b < 0. The angle 8 e appears in a natural way in our 
mathematical analysis and is new in the context of the 
bead and the hoop problem. This angle is represented 
with a star in Figs. |3] and HI 

We call e the "ejection angle" because for any O < 
e , the bead will remain on the initial half of the hoop 
forever and never reach the other side. However, when 
6> > 8 e the bead has sufficient energy to cross the relative 
maximum at = and eject itself from the original half 



of the hoop where it was released. We note that 9 e (oj c ) = 
3 (uj c ) = which compels us to extend their definition 
by imposing e = S = for to < uj c . The two angles 
S and e also satisfy the inequality e > S and e € 
[0, tt] where e (ui — > oo) — > tt. Curiously, we observe that 
V g (0 e ) = 2V g (0 S ) but are unaware if this is the result of 
some physical constraint. Further considerations where 
the bead is placed directly on this point are addressed in 
Sec. lIIIBll 

The integral Eq. ([23]) . along with Eq. (|2"7|) . can be re- 
lated to the first elliptic integral [Eq. . Inversion of the 
integral gives r(i) and thus 0(t) using Eq. Which 
particular inverse JEF equals Eq. (|25[) depends on 71 and 
72 (or a and b). 

For u) < uj c , we have S = e = and a = b > 
for n G [0, tt]. From the tables in Ref. [ID with = a 
and 7I = b in Eq. ([27]) . we identify the cosine-amplitude. 
Transforming back to the variable = 2 arctanr, we ob- 
tain 

-0(ty 



tan 



Ljty/\p-\(a+\b\), 



\b\ 



(31) 



On the other hand, when u > lu c the range 0o G [0, tt] 
is divided into three intervals [0,# s ], [0 s ,9 e ], and [0 e ,7r], 
each of which leads to a solution involving a different 
JEF. For 0o G [0, S ] the roots obey the inequality b > a > 
0; for G [0 s ,0 e ] we have a > b > 0; and for G [0 e ,Tr] 
we have a > > b. Using again the table of integrals 
in Ref. Ojl we set 7^ = a and 7I = b for regions 0q G 
(O,0 S ) and 0o G (0 e ,ir), whereas = b and 7I = a for 
0o G (0 s ,0 e ). This allows identifying the Jacobi elliptic 
functions nd, dn, and cn, respectively. The third interval 
0e < < tt, turns out to lead to an expression that is 
identical to Eq. (|3"Tj) above. 

One can reduce the three cases for uj > lo c to one single 
expression for the entire interval [0, tt] by using proper- 
ties of the JEFs derived in Sec. Ill Al and Appendix [A] 
The simplest form of the solution is found in terms of 
the delta-amplitude function, which naturally appears 
in the second interval. The cosine-amplitude function 
appearing in Eq. (|3ip can be transformed into a delta- 
amplitude with Eq. (|A2|) . Similarly, using the relation 
discussed after Eq. (|A3[) . k' nd(u, k) = dn(u + K,k), it is 
possible to reduce the nd solution in the first interval to 
a delta-amplitude function. Careful consideration of the 
argument (u, k) in the transformation leads to the final 
result 



1 (i) = 2 arctan < y/adn 



-ojty/\p- \ a, 



€[0,tt]. (32) 



This is the exact solution of the equation of motion (|20|) , 
which fully describes the motion of a bead on a spinning 
hoop for arbitrary initial conditions (0o,0o)i excepting 
the cases where i^tot > 1 i / g( 7r ) or w = 0. Equation (|32|) 
simplifies to 0(t) = S for all times if 0q = s . This solu- 
tion improves or complements previous studies of differ- 
ent variants of the problem 1 1 1 2 
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We note that the argument u(t) of the JEF is a linear 
function of time and is proportional to (f> = cut given by 
Eq. (f25j) . As pointed out in Sec. HJ it is only in special 
cases that the argument has a useful physical interpre- 
tation, and the bead on the hoop is not one of them. 
In fact, the problem is slightly more difficult than what 
was presented in the previous section because the JEF 
was only identified after performing a Wcicrstrass sub- 
stitution on the integral of Eq. (|25|) . The Wcicrstrass 
substitution can be shown to be a particular mapping of 
the position of the bead onto the xy-plane. The map is 
obtained from the intersection between the a;y-plane and 
the line passing through the top of the hoop and through 
the bead. The argument u(t) is related to the arc length 
{Ry\P- 1 a /2) 4> 011 this stereographic projection ji£ 



B. Limiting and Special Cases 

We discuss the special cases not considered in the pre- 
vious section: the separatrices (E = 0, 2mgR), the free 
spinning of the bead around the hoop (E > 2mgR), and 
the general solution for the pendulum (u> = 0). For all 
these cases the roots (|2U)) take values a, b —> oo or and 
the determination of 9(t) requires special care. 



1. Separatrices 

Figure [5] shows phase plots for this problem. The sep- 
aratrices form the boundaries between qualitatively dif- 
ferent motions of the system. They can be understood 
as delineating a phase change and typically involve func- 
tions which are asymptotic to some value. Below, we 
obtain simple analytical expressions for the separatrices 
that allow us to plot Fig. [5] without the need of a numer- 
ical procedure applied to Eq. (|25)l . 

As seen in Fig.[5}D, there are two separatrices when u> > 
lo c . The innermost separatrix (dashed thick line) passing 
through (9 0} 9o) = (0,0) separates the cases 9 < 9 e from 
9 e < 6 < 7T. When 9o < #e the periodic motion is 
limited to either of the two wells of V e g (see Fig. 2]) and 
the bead oscillates exclusively on the half of the hoop 
where it started. By contrast, starting positions in the 
range 9 e < 9q < ir involve the bead oscillating between 
the two sides of the hoop and moving through both wells. 

The expression 9(9) of the inner separatrix is obtained 
from Eq. (|2~Tj) by noting that 6q = 9q = is one point 
of the separatrix. At that point C = 1 — 2 cos 9 S , which 
implies that any other point of the separatrix satisfies the 
equation 



= (1 - cos6>) (1 + cos (9 - 2cos6» s ) 



(33) 



•OS 



•D 




FIG. 5. (Color online) Phase plots drawn from analytical 
expressions for (a) subcritical (u> < uj c ) and (b) supercrit- 
ical (uj > toe) cases with the values given in Fig. [4] The 
stars (±0 e ) and diamonds (±0 S ) indicate the stable dynamical 
equilibrium points. The solid lines correspond to the periodic 
motion of the bead inside one of two potential wells. The 
dashed lines are motions of the bead across the two potential 
wells. The dotted lines correspond to the spinning of the bead 
around the hoop. The two thick lines (solid and dashed) are 
the separatrices. 



The motion on the separatrix is obtained from 
Eq. (|2"5)l . The result obtained for arbitrary initial condi- 
tion on the separatrix is rather involved when (9o, 9q) ^ 
(# e ,0) and will not be needed in the following. For the 
simpler case of 9q = 6 e , Eq. (|2l)|) or (|3"2"j) becomes 



The inner separatrix has the shape of a lemniscate as 
shown by the thick dashed line in Fig. [SJd. 



6 {t) = 2 arctan 



%- x scch 



,(34) 
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which is valid for w > oj c . We note that in this case the 
point (6, 9) = (0, 0) is an unstable equilibrium position. If 
the bead starts at any other initial condition (for example 
at 9 e and 9q = 0) the bead will need an infinite amount 
of time to reach the bottom of the hoop. 

The second separatrix divides the phase plot between 
oscillations of the bead that reach both halves of the 
hoop, and the "free" rotation of the bead around the 
hoop when E to t > Kff ax - The equation of this sep- 
aratrix can be obtained from Eq. (|21j) by noting that 
(9q = 7T, 8q = 0) is one point on the separatrix. At that 



point Eq. (f22j) reads C = 1 + 2cos6* s , which immediately 
implies from (|21[) that at any other point of the separatrix 
we have 

s = (1 + cos0)(1-cos0-|-2cos0 s ). (35) 

This separatrix is the thick solid line in Figs.[5ji and[5j). 
Any initial condition away from (9o,0o) = (tt, 0) on the 
separatrix will lead the bead to tend toward 7r for t — > oo . 
To obtain 9(t) for initial conditions lying on this second 
separatrix we integrate Eq. ([25| to get 



6 (t) = 2 arctan 



cos 9 



sinh < u)tyl + cos 9 S + arcsinh 
cos 9 S 



tan - 



' 1 + cos 9 S 
cos 6*c 



(36) 



The results of this subsection leads to the important 
observation that the separatrices are not described in 
terms of JEFs. Instead they depend solely on trigono- 
metric and hyperbolic functions. A further observation 
is that the time evolution along the separatrices is not pe- 
riodic but is asymptotic towards an unstable extremum. 



2. Free spinning 

For initial derivatives 9q greater than the right hand 
side of Eq. (f3"5"l) . we arc outside the second separatrix 
in the phase plot and the bead is spinning freely on the 
hoop. Such motion is obtained, for example, when 9q = tx 
and 9o ^ for which E tot > V^ ax . Any other initial 
condition (9q,9o ^ 0) such that the total energy of the 
bead is larger than V^ ax can be mapped to this initial 
condition at 9q = ir. To obtain the motion 9(t) of the 
bead in this special case it is, once more, easier to start 
with Eqs. ([211)-®. For (#0 = tt, 6 ^ 0) we have from 
Eq. ([26]) p- = 9%/uj 2 , p + = Acos9 s + p- and r -> oo. 
The solution is written in terms of the tangent-amplitude 
function tn = sn/cn. The integration of Eq. (f27f leads to 



9 (t) = 2 arctan 



tn - 



7-<>o 



t. 



71 



but with the roots 71/2 = 7+/- given by 

r 



(37) 



7± 



r 



±2Jr 2 + ^(l-cos6g, (38) 
V 2# 



with r = (w 2 /#o) (1 + cos6> s ). The oscillation of the bead 
corresponds to the open dotted lines outside the second 
separatrix in Fig. [5]) (and outside the separatrix of Fig.[5K 
for ui < LUc). 

For <? — > 00, 7± = — 1 and Eq. (|23|) reduces to the 
form of the arctangent function given in Eq. ((TJ). This 
implies that as 9q increases without bound, our solution 



approaches the form 9(t) 



t. We recognize this as 



characteristic of a free-particle solution. 



3. General solution for pendulum (ui = 0) 



The pendulum is described by Eq. (|20[) for ui = 
and thus reads 9 = — oj 2 sin 9. The solution derived in 
Sec. IIII Al is inapplicable to this case since we divided 
Eq. (f2TJ|) by uj to get Eq. (|2"Tj) . We can nevertheless follow 
a similar path as in Sec. IIII Al for this simplified equation 
of motion. The effective potential (TTS)) with ui = gives 
V c ff = V g {9) and displays a single minimum at 6 = 
(dashed curve in Fig. @|. The equation for 9 (t) with 
O = corresponding to Eq. (|2"Tj) is 



(39) 



9{t) -2uj 2 c cos 9 (t) =C, 



with C = 9q/2 — u! 2 cos 9q. This leads to a phase plot 
similar to that shown in Fig. [5^,. For initial conditions 
located inside the separatrix there are two possible states 
of the bead. For #0 = we have 9{t) = 0, the equilib- 
rium position. For 9q € (0, n) and considering first the 
case #0 = the solution obtained by integrating Eq. (f3"5)) 
(see Sec. IIII A|) involves a delta-amplitude Jacobi elliptic 
function 



9 (t) = 2 arctan < 7+dn 



oj„t cos 



Ih 




(40) 
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with 



1 



COS On 9 
71 



1 + COS On 



(41) 



This expression is a special case of thegeneral solution 
and has been discussed in Refs. 
The scparatrix is obtained from Eq. and (9o = 
with On = ir. Integrating that expression one more time 
leads to the motion of the bead along the separatrix 



(t) — 2 arctan < sinh 



ui c t + arcsinh tan 



(42) 

Finally, the free spinning case can be obtained in ways 
similar to the previous sections. We assume (On = tt, Oo_^ 
0) and obtain Eq. ([57)1 with the new roots (see also Ref.@) 



-1, I 1 - 



1ft 2 

4^0 



(43) 



IV. CONCLUSION 

We introduced Jacobi elliptic functions using integral 
inversion and discussed their properties and the mean- 
ing of their arguments (u, k). The first argument of the 
JEFs, u, was shown to be related to an arc length in 
elliptic geometry, while the second argument k is the ec- 
centricity. As an application of JEFs, we discussed the 
motion of the bead on the hoop and derived the complete 
analytical solution of the problem. Our results reduce 
to the solution of the pendulum when the hoop rota- 
tion frequency vanishes. Interestingly, all motions of the 
bead on the hoop are written in terms of a Jacobi elliptic 
function, except for equilibria positions and scparatri- 
ces, where JEFs reduce to trigonometric and hyperbolic 
functions or constants. The knowledge of JEFs allows 
expressing the complete solution of the bead on the hoop 
problem in an elegant and concise form. These functions 
and their salient features should be remembered as they 
may appear in number of other problems in physics. 



ACKNOWLEDGMENTS 

We gratefully acknowledge the support of the National 
Science Foundation under grant DMR-0907242, the block 
grant of CNSM and the graduate research fellowship at 
CSU Long Beach. 



Appendix A: Additional properties of Jacobi elliptic 
functions 



We introduced Jacobi elliptic functions in Sec. |TT] and 
presented their main properties. We summarize in this 
Appendix a few additional properties of JEFs that were 
used in solving the bead on the hoop problem. 

The JEFs can be related to each other by an appropri- 
ate change of variables. For example, starting with the 
equation u = F(y, k), its inversion gives y = sn(u, k). On 
the other hand, multiplying both sides of the equation by 
k and changing variable t' = kt in F(y, k), we can again 
invert the equation to obtain ky — sn(fcu, fc" 1 ). Thus, 
we found the first of the following relations 



k sn(u, k) = sn(k u, k 1 ), 
dn(u, k) = cn(k u, fc _1 ). 



(Al) 
(A2) 



The second relation is obtained from the first, using 
Eqs. © and (fTU|) . We use the second equation in Sec. lIIII 
Note that in other references, k 2 = m is used. 

The derivatives of the JEFs with respect to u can be 
determined by using their definitions as well. For exam- 
ple, from Eq. (|6]) one obtains dsmi/du = cnudn-u. From 
the Taylor expansions of the trigonometric functions one 
can also expand the JEFs in powers of u. Addition for- 
mulae can be derived from Eqs. (jS])-©. For example, 



sn(u + v) 
dn(u + v) 



snw cm; dnv + sni> cnu dnu 
A(u,v) ' 
dnit dnv — k 2 smt smi cnu cm; 

A(u,v) ' 



(A3) 
(A4) 



with A(u,v) = 1 — fc 2 sn 2 u sn 2 ii. This definition relates 
to A(u) introduced in the paragraph preceding Eq. flSJ 
via the equality A(it) = A(it, K). These addition formu- 
lae allow establishing relations between different JEFs 
as well. We used one of them in Sec. IIII Al namely 
/c'nd(u, k) =dn(u + K, k). This relation can be obtained 
using Eq. (|A4j) with v = K. 

Finally, the JEFs were defined for u G R, but their 
definition can be extended to the complex plane with 
the help of the above addition formulae. For example, 
using Eq. (|A3|) one can decompose sn(a: + iy) into real 
and imaginary parts. In the complex plane, the JEFs are 
double-periodic functions with periods given by integer 
multiples of the complete elliptic integral of the first kind. 

Further definitions, properties, and expressions can be 
found in Refs. S, ElrllE and [3 
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